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04 ' The mean electromotive force that occurs in the framework of mean-field magnetohydrodynamics is studied for 
. cases in which magnetic field fluctuations are not only due to the action of velocity fluctuations on the mean magnetic 

■ field. The possibility of magnetic field fluctuations independent of a mean magnetic fleld, as they may occur as a 
consequence of a small-scale dynamo, is taken into account. Particular attention is payed to the effect of a mean 

Oh rotation of the fluid on the mean electromotive force, although only small rotation rates are considered. Anisotropics 
\ of the turbulence due to gradients of its intensity or its helicity are admitted. The mean magnetic field is considered 
to be weak enough to exclude quenching effects. A r-approximation is used in the equation describing the deviation 
' of the cross-helicity tensor from that for zero mean magnetic field, which applies in the limit of large hydrodynamic 
' Reynolds numbers. 

^ ' For the effects described by the mean electromotive force like a-effect, turbulent diffusion of magnetic fields etc. 
^ in addition to the contributions determined by the velocity fluctuations also those determined by the magnetic field 
J — fiuctuations independent of the mean magnetic field are derived. Several old results are confirmed, partially under 
00 I more general assumptions, and quite a few new ones are given. Provided the kinematic helicity and the current 
C^l ' helicity of the fluctuations have the same signs the a-effect is always diminished by the magnetic fluctuations. In 
0^ \ the absence of rotation these have, however, no influence on the turbulent diffusion. Besides the diamagnetic effect 
. due to a gradient of the intensity of the velocity fluctuations there is a paramagnetic effect due to a gradient of the 
^ii*- ' intensity of the magnetic fluctuations. In the absence of rotation these two effects compensate each other in the case 
-««^ , of equipartition of the kinetic and magnetic energies of the fluctuations of the original turbulence, i.e. that with zero 

■ mean magnetic field, but the rotation makes the situation more complex. The x J-effect works in the same way 
with velocity fiuctuations and magnetic field fluctuations. A contribution to the electromotive force connected with 

Q the symmetric parts of the gradient tensor of the mean magnetic field, which does not occur in the absence of rotation, 
^ \ was found in the case of rotation, resulting from velocity or magnetic fiuctuations. 

^ ^ The implications of the results for the mean electromotive force for mean-field dynamo models are discussed with 
special emphasis to dynamos working without a-effect. 

The results for the coefficients defining the mean electromotive force which are determined by the velocity fluc- 
tuations in the case of vanishing mean motion agree formally with the results obtained in the kinematic approach, 
specified by second-order approximation and high-conductivity limit. However, their range of validity is clearly larger. 
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1. INTRODUCTION 



The mean-field approach proved to be very useful in studying dynamo processes in turbulently moving electrically 
conducting fluids (see, e.g., Moffatt 1978, Parker 1979, Krause & Radler 1980, Zeldovich et al. 1983). The crucial 
point of this approach is the mean electromotive force which is determined by the fluctuations of the fluid velocity 
and the magnetic field. It describes several physical effects like a-effect, turbulent diffusion of the magnetic field or 
turbulent diamagnetism. In many investigations on the kinematic level the magnetic fluctuations are understood as 
caused by the action of the velocity fluctuations on the mean magnetic field. This implies that they vanish if the 
mean magnetic field does so. The mean electromotive force can then be considered as a quantity determined, apart 
from the mean velocity, by the velocity fluctuations and the mean magnetic field. 

There are, however, many realistic cases with magnetic fluctuations which exist independent of the action of the 
velocity fluctuations on the mean magnetic fleld and do not vanish if it does so. We recall here the numerous 
investigations showing the possibility of small-scale dynamos with a zero mean magnetic field (see, e.g., Kazantsev 
1968, Meneguzzi et al. 1981, Zeldovich et al. 1990, Nordlund et al. 1992, Childress & Gilbert 1995, Brandenburg et al. 
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1996, Rogachevskii & Kleeorin 1997, Klccorin ct al. 2002b). As pointed out already by Pouquet et al. (1976) in such 
cases the mean electromotive force has in addition to the contributions mentioned above, which can be ascribed to the 
velocity fluctuations, also others due to those magnetic fluctuations which exist independent of the mean magnetic 
field. 

Several investigations on the mean electromotive force comprising both kinds of contributions have been carried 
out so far. We mention in particular those by Kichatinov 1982 and by Vainshtein and Kichatinov 1983, in which 
an isotropic turbulence with a scale-independent correlation time was assumed. Results concerning the a-effect, the 
turbulent magnetic diffusivity and the turbulent diamagnetism or paramagnetism for an originally isotropic turbulence 
subject to mean rotation have been derived using a modified second-order correlation approximation by Kichatinov 
1991, Riidiger & Kichatinov 1993 and Kichatinov et al. 1994. 

In this paper we present an approach to the mean electromotive force which reproduces such results under more 
general assumptions and reveals new ones, which are of particular importance for astrophysical applications. We 
exclude mean motions of the fluid other than a rigid body rotation with some small rotation rate. We further consider 
only a weak mean magnetic field so that its energy density is small compared to the kinetic energy density. In this 
way we do not consider quenching effects, i.e. reductions of the coefficients defining the mean electromotive force 
with growing mean magnetic field, which are very important in view of the nonlinear behavior of dynamos and have 
been discussed in a number of papers (see, e.g., Gruzinov & Diamond 1994, Cattaneo & Hughes 1996, Scehafcr 1996, 
Kulsrud 1999, Field et al. 1999, Rogachevskii & Kleeorin 2001, Kleeorin et al. 2002a, Blackman & Brandenburg 2002). 
Finally we use a r-approximation in the equations describing the deviation of the cross-helicity tensor from that for 
zero magnetic field. In contrast to the often used second-order correlation approximation it does not totally ignore 
higher than second-order correlations but considers their influence in some summary way. The T-approximation in 
that sense applies in the limit of high hydrodynamic Reynolds numbers. 

In Section 2 we will explain the concept of mean-field magnetohydrodynamics for a homogeneous incompressible 
fiuid. In Section 3 we introduce a Fourier representation of the velocity and magnetic fluctuations, deflne correlation 
tensors and express the mean electromotive force by the cross-helicity tensor. After giving an equation for this 
tensor in Section 4 we introduce in Section 5 the mentioned T-approximation, which leads to closed equations for 
that part of the cross-helicity tensor which determines the mean electromotive force. In Section 6 general relations 
for this mean electromotive force are given, and in Section 7 more specific relations for the case for small rotation 
rates. In Section 8 we restrict ourselves to the limit of weak mean magnetic fields so that the mean electromotive 
force can be considered as linear in the mean magnetic field and specify the correlation tensors for the "original" 
turbulence, that is, the turbulence for zero mean magnetic field. For zero rotation they are determined by simple 
assumptions concerning the deviations from a homogeneous isotropic turbulence, and the influence of a slow rotation 
on the velocity fluctuations is calculated by a perturbation procedure, again based on a r-approximation. We further 
introduce Kolmogorov-type spectra of the relevant quantities. On this basis we deliver results for the coefficients 
defining the mean electromotive force and discuss them with special attention to the different effects of velocity and 
magnetic field fiuctuations. Finally, in Section 9 our results are compared with results of the kinematic approach 
in the second-order correlation approximation, some remarks concerning their range of validity are made, and some 
prospects are mentioned concerning the extension of the approach of this paper to related questions of mean-field 
magnetohydrodynamics. 



2. FORMULATION OF THE PROBLEM 



We consider a turbulent motion of an electrically conducting incompressible fluid in interaction with a magnetic 
field. Let us assume that the fluid velocity U and the magnetic field B are governed by the equations 

aU , , VP 1 , 

— + U • V U = + — V X B X B 

at p fip 

+ z^AU + 2U X rj + F , (1) 
— = V X (U X B) + r?AB , (2) 
V-U = V-B = 0, (3) 

where P is a modified pressure and F an external force. We refer to a rotating frame with O being the angular 
velocity responsible for the Coriolis force. As usual i' and rj are the kinematic and magnetic viscosities, p is the mass 
density and p, is the magnetic permeability of the fluid, all considered as constants. 

Wc further assume that there is an averaging procedure which defines for each quantity Q an average (Q) and 
satisfies the Reynolds rules. In the spirit of the mean-field concept we split the fiuid velocity U and the magnetic 
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field B according to 

U = U + u,B = B + b, (4) 

into mean parts, U = (U) and B = (B), and fluctuations, u and b. Analogously we split P into P and p, and F into 
F and f , etc. The mean fields U and B satisfy the equations 

^ + (U-V)U = -— + — (V xB) xB + i/AU + 2Ux f2 + :F + F, (5) 

ot p pp 

— = V X (UxB + £)+7?AB, (6) 
at 

V-U = V-B = 0, (7) 

with ^ and £ being a pondermotive and an electromotive force due to fluctuations, 

:F = -((u- V)u) + — ((V xb) xb) , (8) 

pp 

£ = (u X b) . (9) 
For the determination of !F and £ we need information on the fluctuations u and b. These have to obey the equations 

^ = -(U-V)u-(u-V)U-^ + — (VxB)xb 

at p pp 

+ (V X b) X B] +T + i^Au + 2u X n + f , (10) 

— = V X (u X B + U X b) + 77Ab + G , (11) 
at 

where T and G summarize terms nonlinear in u and b, 

T = ((u • V)u) - (u • V)u + — [(b X (V X b)) - b X (V X b)] , (12) 

pp 

G = V X (u X b - (u X b)) . (13) 

3. FOURIER REPRESENTATION, CORRELATION TENSORS 

Let us now represent quantities like u and b, or U and B, by Fourier integrals defined according to 

Q(x,i) = y"Q(k,i)exp(ik-x)d^fc . (14) 

We may rewrite the equations (|l^) and ( pi] ) for u and b into equations for their components Ui and bi with respect to 
a Cartesian coordinate system, derive equations for their Fourier transforms Ui and hi and subject the equation for Ui 
to a projection operator Pij = Sij — kikj/k'^, with Sij being the Kronecker tensor, in order to eliminate the pressure 
p. In this way we obtain 

du,ik)/dt = -Py[4(u,U;k) + S'j(u,U;k)] + (2P„- - S^j)Lj{h,B;k)/ pp 

+S,{h,B; k)/pp + A,(k)u,(k) - j^k^Mk) - P.,(f,(k) - fj{k))/p , (15) 
dkik, t)/dt = 5,(u,B;k) -5.(b,U;k) -L,(u,B;k) + L,(b,U;k) 

^r]k%{k) + G,{k) , (16) 



where 



L,(a,A;k) = i J aj{k ^ K)KjA^{K) (fK, (17) 
5',(a,A;k) = ik^ [ a^{k ^ K)Aj{K) d'^K , (18) 

u2 



A,(k) = 2£,,„fc„(k.^2)/fc^ (19) 
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and Sijk is the Levi-Civita tensor. For the sake of simphcity the argument t is dropped everywhere. Concerning the 



derivation of (|19| ) we refer to (Al) 



We wiU use these equations for calculating the two-point correlation functions. Let us consider, e.g., the correlation 
tensor (wi(xi)?«j (X2)) of two vector fields v and w, where xi and X2 denote two points in space but both fields are 
taken at the same time. Using the definition (Q) of the Fourier transformation and following a pattern introduced 
by Roberts and Soward 1975 we write 



where 



(wi(xi)u;j(x2)} = J y (u,;(ki)wj(k2)) expi(ki • xi + k2 • X2) d^fci d^fc2 
^ J (p,j (r, K) exp {iK ■ R) dPK 

= j ip,j (k, R) exp (ik ■r)d^k , (20) 

(^,j(r,K) = J{v,{k + K/2)wj{-k + K/2))exp{ik-r)(fk, (21) 

(/3y(k,R) = J{v,{k + K/2)wj{-k + K/2))exp{{K-R)(fK, (22) 

and R = (xi + X2)/2, r = xi — X2, K = ki + k2, k = (ki — k2)/2. We relate later r and k to small scales and 
R and K to large scales in the physical space. 

In the following we consider in particular the correlation tensors for the velocity and magnetic fluctuations, 
(Mi(xi)uj(x2)) and (foi(xi)6j(x2)), and the cross-helicity tensor, (ui(xi)6j(x2)), and we use the definitions 

Wy(k,R) = $('Uj,'Uj;k,R) , (k, R) = $(6,, 6^; k, R)/^p , (23) 

X„ (k,R) = $(w„S,;k,R) , (24) 

where 

$(?;, w; k, R) = J {v(k + K/2)u;(-k + K/2)) exp (iK ■ R) d^K . (25) 
The definition of Vij implies 

Wy(k,R) =«,,(-k,R) (26) 

and 

v,j{k, R)k, = ^y^v^J (k, R) , (k, R)k^ ^ Vjf„- (k, R) . (27) 



Here and in the following V stands for d/dR. A relation analogous to ( p6[ ) applies to , too, and relations analogous 
to ( PtI ) apply to rriij and even to Xij ■ 

If we know Vij, rriij and Xij we may calculate the pondermotive force ^ and the electromotive force £. In this 
paper, however, we will focus attention on £ only. Since 

E,=e^jk j Xjfe(k,R)d3fc (28) 



we first deal with the cross-correlation tensor 



Xij- 



4. THE EQUATION FOR THE CROSS-HELICITY TENSOR 

According to the definition ( p5| ) of Xij its time derivative is given by 

'-M^ = R) + Hu.. §; k, R) . (29) 
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We evaluate this using equations (15) and (|6|) and neglecting all contributions containing higher than first order 
terms in the operator V, independent on whether it acts on w^, to^ , Xiji U or B. With manipulations explained in 
Appendix B we obtain 



dt 



L. 



ijklm 



+ Mijklrn^ mXkl + NijklXkl + Ci- 



(30) 



where 



ijklm 

Nijki 



^i{Dik^jlkm ~ £ifcm'5j7(k • O) — EikpkpSjli^m) / k 

+i{ri - iy)SikSjikm , 

U^^kSji - SikUjA ~ 2kipUp^kSji - DikSji + {v + ri)k'^StkSji 
-i(k ■ B){v,j - niij) + -;^BpSIp{vij + m^j) 



(31) 

(32) 
(33) 
(34) 



(35) 



and Bi,j — V jBi, Ui.j — S/jUi, furthermore kij = kikj/k^, Vijq = ^dvij/dkq and niijq = ^drriij/dkq. 



5. THE r-APPROXIMATION 

When dealing with equation ( pO| ) for Xij we are confronted with the difficulty that we know nothing about the term 
Cij which stands for couplings of Xij with higher-order correlation tensors. We assume now that we know the solution 
Xij of equation ( |30| ) for vanishing mean magnetic field, that is lij = 0, denote it by Xij\ ^'H'i the corresponding dj 

by C^j"^ . For the case with non- vanishing mean magnetic field we put then Xij = Xif' + Xif' ^-nd Cij = C^^"^ + c\^\ 
Then we have 

, T (^^kl , V7 (B) (B) ^(B) _ J ,„„^ 

+iv,jfcim-g^ V M^jkhyiy mXkl +'^^OklXkl + " • (-^OJ 

At this level we introduce the r-approximation 

c{f^=xlfMfc), (37) 

with some relaxation time r(fc) (see, e.g., Orszag 1970, Monin & Yaglom 1975, Pouquet et al. 1976, McComb 1990, 
Kleeorin et al. 1990, 1996). This approximation applies for large hydrodynamic Reynolds numbers, i.e. for fully 
developed turbulence. In this case the relaxation time T(k) is determined by the correlation time of the turbulent 
velocity field only. The limit r — > 00 corresponds to cancelling all higher-order correlations, that is, to some kind of 
second-order correlation approximation (see also Section 9). Note that we do not introduce any closure assumption 
for the original turbulence; equation (^) concerns only deviations from the original turbulence. 



6. RELATIONS FOR THE MEAN ELECTROMOTIVE FORCE 

Returning now to the electromotive force £ we first put in the above sense £ — + £''^\ We assume however 
that £''^\ if not at all equal to zero, becomes unimportant in comparison to when B grows and put therefore 
£(0) _ Thus the electromotive force £ is given by 

£^ = e,,k J xj-f (k,R)d'fc. (38) 

Let us assume that the time variations of lij, that is, of Bi, Bij, Vij, and rriij, are sufficiently weak, so that xij 7 
apart from some initial time, can be considered as a solution of the steady version of equation (^6|), that is, with 
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dXij /dt being neglected. Thus the dependence of xij on Bi and Bi j is hnear and homogeneous and, in addition, 
instantaneous. So we may conclude that 

= o^ijBj + bijkBj^k (39) 

with tensors Gij and bijk determined by Vij and rriij. 

This relation can be rewritten in another form, which might be more suitable for some discussions. To derive it we 
split Uij in a symmetric and an antisymmetric part and express the latter by a vector. We further split the gradient 
tensor of B into a symmetric part c?B, defined by [dB)ij ~ ■^{Bi.j + Bj ^), and an antisymmetric one, which can 
be expressed by VxB. We have then Bij = {dB)ij — ieijfe(VxB)fe. Finally we express the tensorial coefficient 
occurring then with VxB again by a symmetric tensor and a vector. In this way we arrive at 

£ = -aB-7xB -/3(VxB) - Jx(VxB) ndB (40) 

(see Radler 1980, 1983), where a and (3 are symmetric tensors of the second rank, 7 and S vectors, and k is a tensor 
of the third rank; k can be considered to be symmetric in the indices connecting it with c)B, and contributions can 
be dropped which would produce V • B. We have 



\^{aij + ciji) , Pij — ^( 

1 . 1.. , s 1 



c^ij — ^77(0^" + 'iji) I Pij — -ri^ikibjki + £jkibiki) , (41) 



7i — 2^ijkO"jk , — ^{^jji ~ bjij) 7 ^^ij k — — -^{bijk + bikj) ■ (42) 

For the sake of simplicity we restrict ourselves on the case Ui — const. Thus xlj is governed by the steady version 

of equation (^) with Lijkim = and Nijki = -D^kSji + + ifjk'^SikSji- Using (|^ the equation for Xif ^ can be 
written in the form 

DikXkj + Mijklm'^mXkl = ^v'^* ' (43) 

where 

Aj = - D,jT^ , ^ + + T])k^ . (44) 

In deriving equation ( p30| ) we have neglected all contributions to lij of higher than first order in V. In the same sense 
we may replace '^mXa equations (^0|), ( ^6| ) and (43) by VmXfci^\ where Xij^^ summarizes the contributions to 
Xij which are of zero order in V. We have then 

D^kxlf = , /° = -*(B • k)in,, - m.,) , (45) 

and consequently 

xT - At'^S^* , (46) 
where is the inverse of D, satisfying B^^Bkj — that is, 



U2 



A/ = (1 + t^^fc^)"^(% + t^e»jfefcfc + u?hkj) , (47) 



where u: — 2T*(r2 • k)/fc . Proceeding as described above we finally obtain the solution of equation (|4^) in the form 

- DT-^'hjn - D^p'M„klmD^,'VraI°lT^ • (48) 

7. THE MEAN ELECTROMOTIVE FORCE AT SLOW ROTATION 

We split the electromotive force £^ ' according to 
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into parts f^^"^ and £^^^\ the first of which does not depend on the rotation while the second one does but vanishes 
with vanishing rotation. For e[^'^^ we have simply 

Sf'^ = e^k f I 3k (k)T, (fc) (fk . (50) 



As for we restrict ourselves to the case of slow rotation. More precisely, we neglect terms of third and higher 

order in fir* in comparison to unity. So we find 



kik- 

'4(k • nfie^jk + em-^)IjkT* 



-2z(k . n)(2(i^^(k . V) - (O . V))e.,./]:) 

-^^£.,fefc,V,4°V*)§d'fc. (51) 
In the integrand we have dropped terms of the structures Vjfcfe/^°^ or Vjfcfc//^\ which would, by reason connected 



with (^7|), lead to contributions of the second order in V. Moreover, for the sake of simplicity in both (^0|) and (|51^ 
we ignored terms containing the factor (77 — i^)fc^T». We will later point out the consequences of that for our final 
result. 

In relation ( p9| ) for £ we split a.y and into parts corresponding to their dependence on vij or m^j, e.g., 

a.,=a!^+o^. (52) 
and each of them are splitted into one which is independent of and a remaining one, e.g., 

«r-«r+«r^- (53) 

We note that the splitting of S^^'^ according to ( ^9| ) and the splittings of a|j\ a|™', • • • according to ( p3| ) refer only 
to the dependencies on Q, which occur explicitly in this stage of our derivations. In general will depend on Jl, too, 

and then 5^^") and a|j°\ 0'^ij^'^\ ' ' ' '^ill have also contributions with $1, and the dependencies of S^^^^^ and a-J^\ 
^(msi) ^ . . . j-j ^jjj more complex. 
A straightforward calculation yields 

a-J°^ = eii^v\2j , (54) 

= -e^l„,mhlj , (55) 

C - ^'.'-/^^ ' (56) 

b^^k^ = eyi"^lfe^ + 2e»im"iS„ , (57) 

-4£.p««S/™f^i^^m , (58) 
^(mn) ^ (-2mJ,p',^., + V,mJ,J^., + VjmJ,^, - 2Vp(to21^-,p + mg];) + V/mJ,^^.)^; 

(3) 



+4£vq"l),^j7m^i^m , (59) 
"jjA; — K'^^ikjl ~ '^^jkil + Vppijkl>^''l 

-^e^jpi'L^i^^ , (60) 

-4(£ypmg„^ + 2e^pgm^J^.,,„Jf],f]™ . (61) 
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Here we used the definitions 

vifk^.p ^ y'(-*fc)%(k)fc,...pr^'d3fc, (62) 

vljklmn = J V^.jik)kkhnnidT^/dk)k(fk, (63) 

where A = or A = 1 for even or odd v, respectively, and kk...p = kk---kp/k'^ , with v being the rank of this last 
tensor, and analogous definitions for fn'^jl p, and m'^jkimn with Vij replaced by niij. With the help of ( p6| ) we find 
that the v^^l ^ with A = are symmetric and those with A = 1 antisymmetric in i and j and moreover in both 
cases symmetric in every pair of the remaining indices, and that v'^jj^i^^ is symmetric in i and j and again in in every 

pair of the remaining indices. The same applies to m'^^jl ^ and to rn'^-j^i^^. Furthermore in the case A = we may 

conclude from (|7|) that w^y- „p = ^(l/2)VpUp^j „ and w-pj „p = (V2)VpU^pj „, where the p (which will never 

explicitly occur in the following) are defined like the v\'^^ ^ but with kij,,,pk~'^ in the integrand instead of kij p. In 

the case A = 1 we have simply Wp^] „p = (l/2)VpUp^j „ and w-pj „p = -(l/2)VpW,[pj „. AU these statements apply 

again to Wp^j „p and m^pj^^^^p, too. We have used these properties of the v'^^l ^, '^[fk...p^ deriving (|^~(|l]), 

and we have ignored all terms which would result in contributions to £ of higher than first order in V. Likewise we 
have cancelled contributions to biji^ proportional to djk^ which can not contribute to £ since V • B = 0. 



Let us consider the representation (40) for £ and split each of the quantities a, 7, f3, S and k after the pattern of 
(b^) and (|5^) into four parts, e.g.. 



a!;^ + , (64) 

(f) _ (fO) (vCl) (m) _ (mO) (mQ) 



and 

= "y- + alj , a-j = + a]^ . (65) 



The above remarks concerning the dependencies on fl apply here analogously. Using (p])-(p2[) and (^4|)-(^l]) all these 

contributions to a, 7, /3, S and k can be expressed by the v^^l p, m[jl p, • • ■. For the sake of simplicity we give here 
only the contributions which do not depend on fl. They read 

(66) 
(67) 
(68) 
(69) 
(70) 



(i;0) 
Q;- ■ — 


^1 (1) 1 


- £jifc4l]) , = -\i^akm^kl + 






2^ vv y 




- 2-i;u , 




2 ^1 ' 




= , 


JvO) _ 
1^1] k - 


-^^^^ji^ik 






(mO) _ 




+ eikim[]^) + Siimim'^kiji + mj^fc,) 





In order to derive (|67|) we have used (A2) 



8. THE MEAN ELECTROMOTIVE FORCE IN THE LIMIT OF WEAK MEAN MAGNETIC FIELDS 
UNDER SPECIFIC ASSUMPTIONS ON THE ORIGINAL TURBULENCE 

In general the correlation tensors Vij and rriij depend, of course, on the mean magnetic field B. As a consequence 
the mean electromotive force £ depends in a nonlinear way on B. We restrict ourselves now to the approximation in 
which Vij and rriy are simply replaced by the corresponding tensors v^^^ and for the "original" turbulence, that 
is the turbulence for zero mean magnetic field. Then, of course, £ is linear in B. This implies that any quenching 
effects are excluded. 
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8.1. Nonrotating turbulence 



We first ignore any influence of a rotation of tlie fluid on the turbulence. As for the correlation tensors v^^"^ and 
m^ij"^ for the original turbulence we assume that, as long it is homogeneous, they are essentially determined by the 

kinetic and magnetic energy densities, that is by (u*^°)^) and (b'^*'^^), and by the kinematic and current helicities, 
(u^"^ • (V X u'"')) and (b'"' • (V x b*^"))), and that inhomogeneities are only due to gradients of these four quantities. 
With the notations u^^^ and b*^"^ instead of u and b we want to stress that we are dealing with the original turbulence. 
Note that b is in general non-zero even if b''^^ = 0. Under the assumptions adopted the most general form of the 
correlation tensor v^^^ (k, R) of the velocity fluctuations of the original turbulence is given by 

-^[e^jkkki2i + ^(k • V)) - ^{hejim + kjeiim)ki\7^]^i^''\k,-R)^ . (71) 

We note that satisfies the requirements resulting from V • u^") — 0, and that VF^"'(fc,R) and ^^''^(A;,R) 
are spectrum functions depending on k via k only and possessing the properties W^'"\k,'R) dk — (u^''-'^) and 
fi'^'"^ (k^H) dk = (u'"' ■ (V X u'"')). The definition of the correlation tensor m^°^ of the magnetic fluctuations of 
the original turbulence follows from ( [7l| ) when replacing W'''"^ and fi'^'"^ by mf^\ M^^™) and and we have 

/o°"ty(™)(/t,R) dk = (b(o)2)/^p and /g" M(™)(fc, R) dfc = (b^ • (V x b(o)))//ip. 



Since we ignore here any influence of rotation the a^J , ct^™ , • • • coincide with the , a^™ , • • • given with 



(|66|)"([70|). Calculating now the w-j^^...p with replaced by vlj as specified by ( [h] ) we find 



i<5,,/(^'i) + • • • , (72) 



vSl - ^(S^kV, 6,,V.,)&^^ - \e,,uJ^^'^ , (73) 



ijkl g V "iJ ft^t ^ ijkl 

where 



= iy('')(A:,R)T,dfc , = / ^(")(fc,R)T,(ifc 

Jo Jo 



(75) 



and E^jj,^ is a completely symmetric tensor, which we define with a view to the following in a more general frame by 

The terms in ( fz^ ) and ( [T^ ) indicated by ■ • ■ are without interest because they would lead to contributions to £ which 
are of second order in V. For the , ^-iid '^ijL analogous relations apply with analogously defined J^™^) and 

j(ml)^ 

Using now (^-([70|) we obtain 

= \5^rJ^^'^ , = , (77) 

= Is^^li^^) , = , (78) 

= iv/('"i) , 7(") = _lv/("i) , (79) 
Siv) ^ = , k(^) = . (80) 
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8.2. Rotating Turbulence 

In order to include the effect of a rotation of the fluid we have first to study how the original turbulence, that is w^^'' 
changes with the rotation, that is with $1. We assume here that in the absence of rotation the turbulence possesses 



no helicity. More precisely, we assume that ( |7l| ) applies with /i^^^ = in the limit of vanishing Jl, and we will now 
calculate the effect of rotation by a perturbation procedure up to the second order in ft. 
Proceeding as in the calculation of the cross-helicity Xij we start from 

9wi,(k, R) ^ ,dui _ , duj , ,„ , 

dt =^(^^"j;i^^R) + ^("'^^;k,R) (si) 

but use the equation ([l5| ) for Ui with B = and T = so that Vij turns into vf^\ and put again U = 0. So we arrive 
at 

^ MyfeimV„fi°' + N.jkiv^^^ + 2vevf''^ + C^, = , (82) 



dt 



where 



Mijklm = ■^^p{PmpO<)kq + Pmq{k)kp){eikq5jl ~ EjlqSik) , (83) 

Ntjki = -2 — p — {eikpSji + ejipSik))kp , (84) 

= -'i>{u,,P,kfj)-<^{Akfk,Uj) ■ (85) 

We assume now that we know the solution v^j'' of (^) for = 0, denote it by w^^^^ and the corresponding Cij by 
ci^\ For the case with non- vanishing Jl we put then wf^-* = uf^"'' + and Ca = (71° + cl^\ So we find 

ij ^ ^ ''J '■J '■J '■J '■J 

— \- M^jklmy-rnVki + '■^ijMVki + Zvk V^j + O,^- -7-^. , («bj 



where 



At this level we introduce again a r-approximation. 



4f ^ = ~{M,,klm^rn + N,,kl)v^^°'^ ■ (87) 



cf) = .fV^(k,i^), 

where f(k, fi)) is a relaxation time analogous to r(fc) introduced with (^). Now we assume that the characteristic 
time of variations of vf^^^ is much larger than f so that we can drop the time derivative dvf^^^ /dt in (86). Thinking 
of slow rotation, that is fir <C 1, wc expand vff^^ in the form 



_ ^,(01) , ^,(02) ^gg^ 



where v^j and v^j are of the first and second order in O, respectively. We further assume that f does not deviate 
markedly from r, that it does not depend on the sign of and therefore its expansion with respect to $7 possesses 
no linear term. In that sense we put f = r + 0(f7^), where the last term is without interest for the following. In this 
way we we obtain 



V. 



= -n(M,,fei™V„, + Nrjki)vfi'' '\ K'=T-'+ 2yk' , = 1, 2 . (90) 



Identifying now w^;"'' with u^^^* defined by (|7^) with /i'") = and neglecting again terms of higher than first order in 
V we find 

- S4^^."fc»(^(k • V) - (O . V)) W^(")(fc,R) , ) = . (91) 
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(0).. 1 „ . i 



With we arrive at 

= 8^(^^^^-« + — (fc.V,-fc,V. 

+^e^„nk^{^^^ik■V)-i^l■V))^W^^Hk,■R) . (92) 

Interestingly enough this coincides with ( [7l| ) if we replace there by /i'"-*, defined by 

il^'"^ = -2f,(^^^j^{k-V)-in-V)^W^-^ik,R) , (93) 

and cancel the terms with V/i'"-', which are of second order in V. Note that the absence of terms of higher than 
second order in fl in the correlation tensor (^2|) is a consequence of the fact that we restricted ourselves to small Jlf 
and used in that sense ( |89| ) . The absence of terms of the second order results from the assumption on the dependence 
of f on fl. 

As for the magnetic fluctuations of the original turbulence we assume that they are independent of the rotation of 
the fluid and show no helicity so that 

mf{k,R) ^^(^P.^ik) + ^(fc,V. - fc.V,))w^(")(fc,R) . (94) 

We calculate now the vljl,,,^, m'ij\...p, ■ ■ ■ defined by (^)-(^) with Vij replaced by according to (^) and find 

^ = + - > (95) 

«S = ^('5»^V,-<5,,V,)/('''^) 

-\e^,,{5A^ ■ V/('''^)) - ^<,Lf^;V™/('''^)) , (96) 

«S = \(^^M~\e^i)I'^'' + -- - , (97) 

-\^^Jp{E%^{^ ■ V/('^'')) - y<L,.f^,V./("^)) , (98) 

>^ - J_fA..s;(4) „:^c;(6) 
In accordance with (|75[) we have used here the definitions 







iM = M/('')(fc,R)T,^(fc)dfc , 7("^) = / w'^''\k,R)Ti:{k)n{k)dk . (100) 



The expression for w-j'^Zmp turns into v'^^j^i^^ if /('"'^) is replaced by /("^' where 

/"OC 

/(^)'= / W^''\k.,-R.){dTl/dk)kdk . (101) 







apply for m\^\ m\^\, ■ ■ •, too, but with /(™'^) = 0. 



Again terms indicated by • • ■ correspond to contributions of second order in V to £. Relations analogous to (p5[)-(101) 



Calculating now the a|j\ " ■ ' according to (|5^)-(|6l|) and finally the a^j\ cx[j'\ ■ ■ • according to (p])-(|4^) we 
arrive at 
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7 

(m 



_ 1 

^10 V3 
2 

6 5 6 15 

6 5 6 15 

6 

6 



6 

30 



(102) 

(103) 
(104) 
(105) 
(106) 
(107) 
(108) 
(109) 
(110) 
(111) 



8.3. Specification to Kolmogorov Type Turbulence Spectra 

Let us now specify the original turbulence to be of Kolmogorov type, i.e., to possess a constant energy flux through 
the spectrum, and consider the inertial range of wave numbers, ko < k < kd, where k^^ = Iq defines the largest length 
scale and fc^^ the dissipative scale of the turbulence. In this range we have = (g - l){(u'^"^^)/ko){k/ko)-'' , and 
W'^"^^ analogously. Furthermore we put /^^^^ = (g - l)((u(") ■ (V x u(")))/fco)(fc/fco)"'^ but, by reasons connected with 
the conservation of the magnetic helicity in the high-conductivity limit, /u*^™-' = {h^^^^ ■ (V x h'^^'>))S{k — fcg) . Finally 
we assume that = f* = 2ro(fc/A:o)^^^, with 2to being a correlation (or turnover) time for k = ko. In a ll case s q is 
a constant constrained by 1 < q < 3. Assuming that it is sufficient to take the integrals in ( [75| ) and ( |lO0| )-( p^ ) over 
the inertial range only and that k^/kd 1 we obtain 



/(■"A') 



2^ 
Ai+ 1 



j(m)/ 



j{v)l 

2^ (bW^)_ 
-2{q - l)/("2) 



-2(g - l)/('^2) ^ 



' 



j(ml) 



(b(") ■ (V X b (°))) 
HP 



-To , 



(112) 
(113) 
(114) 
(115) 



Note that the J^^^), /("i), and are independent of g. 

By the way, with the above specification of we have simply (b'") • (V x b^"))) = (a^''^ •b("))//Q. Here (a(°) •b^'^^} 
is the magnetic helicity, where a^"^ is the vector potential of b^°\ i. e., h^^^^ = V x a^°). The factor 6{k — ko) in the 
function Ai^'') is chosen in order to meet the realisability condition for the magnetic helicity (see, e.g., Moffatt, 1978; 
Zcldovich et al. 1983). 



8.4. Specific Results for Nonrotating Turbulence 



Let us now summarize and discuss our results. When speaking in this context simply of contributions of velocity 
and magnetic field fluctuations we refer always to the origin al tu rb ulen ce and so to u^^^^ and b^°^ Starting with the 
case of nonrotating turbulence we specify now (|7^)-(^) by (112)-(115) and obtain 



= V(u(") • (V X uW)) 



(bW ■ (V X b("))) 
lip 



To , 



(116) 
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6 V MP / 

^ = , K = Q . 

We have an isotropic a-effect, and a is a sum a*^"-' + a*-™-' of two contributions determined by the kinematic heUcity 
and the current hehcity of the original turbulence. Whereas the signs of a^"^ and the kinematic helicity coincide, 
those of a'-™^ and the current helicity are opposite to each other. This is in agreement with results, e.g., by Pouquet 
et al. 1976, by Zeldovich et al. 1983 or by Vainshtein and Kichatinov 1983. The contributions of the kinematic and 
the current helicities to the a-cffect, acting in opposite directions, may even compensate each other. 

We also have an isotropic mean-field diffusivity (3, and this is determined only by the intensity of the velocity 
fluctuations of the original turbulence. There is no contribution of the magnetic fluctuations of the original turbulence. 
This again agrees with results by Vainshtein and Kichatinov 1983. 

An inhomogeneity of the original turbulence leads to a transport of mean magnetic flux which corresponds to that 
by a mean velocity —7. The expulsion of flux from regions with high intensity of the velocity fluctuations has been 
sometimes discussed as "turbulent diamagnetism" ; see, e.g., Zeldovich 1956, Radler 1970 or Krause and Radler 1980. 
Our result shows, again in agreement with Vainshtein and Kichatinov 1983, that magnetic fluctuations act in the 
opposite sense, that is, in the sense of a "turbulent paramagnetism" . The magnetic flux is expelled from regions 
where (u'^*'^^} — (l//xp)(b(°^^) is higher, and pushed into regions where it is lower compared to the surroundings. With 
equipartition of kinetic and magnetic energy, (u'^°)^) = (l//ip)(b("^^), this effect vanishes. 



(117) 

(118) 
(119) 



8.5. Specific Results for Rotating Turbulence 



In the case of the rotating turbulence in which helicity occurs only due to inhomogeneity and rotation we have 



15 V 3 ^^p 

-il(f^*V,+f^,*V.)((u(°)^) + Ai^) 



1 lh^">^) \ 
+ -(e,z„f]ff]* + £,,„f]ff7*)V„((u(°)2) _ A L) U.„ . (120) 

A jip J 



ik-((u(°)^)(l-i^f7*^) + lf^ 



5 5 fip 

-^f]*17*((u(°)2)+3^^)Vo. (121) 



7 = ^ v((u^"^^) - ^ ^)(i - -n*') + -n* X v((u("^^) , 

\ /ip 5 3 pp 

+ ■ V((u(")2) _ S 1)) ro . (122) 

5 fJ-P J 

S = -i((u(o)2)-i^)OVo, (123) 
9 ^ f-ip 

+ %,,intni+em^t^*){{u^°^') - ^ ^) TO . (124) 

pp J 

Here Vt* stands for JItq . Our approximations are justified for \ ^ 1 only. We recall that we have ignored all terms 
containing factors (77 — v)k'^T^ in (|5^) and (^l]). The only influence of these terms on the above result would consist 
in the occurrence of contributions to cty which, in comparison to others, are smaller by a factor of the order o f the 
small quantity (77 — v)kQTQ. As far as contributions of the velocity fluctuations are concerned the result (120)-(124) 
is in qualitative agreement with earlier results (e.g., Radler 1980). 
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The a-effect occurs now as a consequence of the simultaneous presence of a rotation of the fluid and gradients in 
the intensities of the velocity or magnetic fluctuations, and it is clearly anisotropic. A rough measure of the a-effect 
is the trace of aij. For this quantity we have 



_ in*.V)l{u^">')--^ '-TO. (125) 

9 V ^ W J 

At least with this quantity the effect of the velocity fluctuations is again diminished by the magnetic fluctuations. 

Like the a-effect the mean-field diffusivity, too, is in general anisotropic. Interestingly enough, in contrast to the 
case of nonrotating turbulence is now no longer independent of the magnetic fluctuations. For equipartition of 
kinetic and magnetic energy, (u*^"^^) = (l//ip)(b(°)^), the mean-field diffusivity is again isotropic. Even in the general 
case the tensor Pij has no other non-zero elements than diagonal ones. If Jl* is parallel to the xa-axis we have 

3 ' lb ixp 

/333 - i(«(")^)(l-^^^*^)-^^^^*^ (126) 

6 o lb fip 

In the absence of magnetic fluctuations not caused by the mean magnetic field, i.e. b*-"-' = 0, all elements of j3ij 
decrease with growing |r2*|. 

Compared to the case of nonrotating turbulence, the vector 7 describing the transport of mean magnetic flux has an 
additional term perpendicular to the rotation axis. This term contains no longer the gradient of (u'"'^) — (l//i/3)(b("^^) 
but that of (u(0)2) + {1/ pp){h'-0)^), that is, it does not vanish with (u^o)^) ^ {l/pp){h<-°^^). 

In contrast to the case of nonrotating turbulence even in the absence of gradients of the turbulence intensities, the 
vector S is no longer equal to zero. So we rediscover the contribution to the mean electromotive force proportional to 
r2x(VxB) , which has been sometimes discussed as "J7 x J-effect" . We note that x (V x B) = — (fJ- V)B + V(r2-B) 
and that the term V(r2 • B) plays no part in the induction equation for B as long as the coefficient connecting S 
and fl does not depend on space-coordinates. We recall that an electromotive force proportional to f2 x (V x B) 
is, even in the absence of an a-effect, in combination with a differential rotation capable of dynamo action; see, e.g., 
Radler 1969, 1970, 1980, 1986, Roberts 1972 and Moffatt and Proctor 1982. With respect to 6 velocity and magnetic 
fluctuations act again in the opposite sense, and d vanishes with (u^'^^^) — {l/pp){h'-^'>^). For a rotating turbulence, 
again even in the absence of gradients of the turbulence intensities, also K^-fe is unequal to zero. 



We recall that we determined the correlation tensor v^^^ used for the calculation of (120) - (124) under the assump- 
tion that the expansion of the relaxation time f with respect to SI has no linear term. As a consequence contains 
no terms o f sec ond ord er in fJ. It can be followed up easily, however, that a deviation from this assumption would 



change in (|120| ) - ( [124D nothing else than numerical factors of the terms of second order in fl* , and it would leave 



(125) unchanged 



Our results (120) - (|122| ) are different from those given in papers by Kichatinov (1991), Riidiger & Kichatinov 
(1993) and Kichatinov et al. (1994). In the derivations of these authors more restrictive assumptions have been 
used. In particular, they dropped the nonlinear terms in the equations for the velocity and magnetic field fiuctuations 
from the very beginning, which is only valid either for small hydrodynamic and magnetic Reynolds numbers or for 
high magnetic Reynolds and small Strouhal numbers. They further assumed that v — r] = I'^/tc, where v and 77 are 
kinematic viscosity and magnetic diffusivity, and l^. and correlation length and time of the turbulent velocity field. 
This in turn applies only if the hydrodynamic and magnetic Reynolds numbers are of the order of unity. Our results, 
which are derived on another way, apply for large hydrodynamic and magnetic Reynolds numbers. 



8.6. Implications for Mean— Field Dynamo Models 



Let us add some remarks on the possibilities of dynamo action of the induction effects described by (|120| ) - ( 124 ). 
For the sake of simplicity we consider an axisymmetric dynamo model. We use corresponding cylindrical co-ordinates 
r, Lp and z. The mean motion is assumed to consist in a differential rotation, 

U = we2xr, (127) 



where the angular velocity lo may depend on r and z, is the unit vector in z-direction and r the radius vector. 
Further, again for simplicity, only the contributions to £ which are linear in il* and independent of b*^*^^ are taken 
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into account. We write 

£ = -ao(g-r2*)B + ai((0*-B)g+(g-B)rj*) -/3oVxB 

-7ofifxB-7i((n*.B)g-(g.B)f2*) (128) 
■ V)B - V(n* • B)) + ■ V)B + V(n* • B)) 

where ao - ffai = if (u(o)2)ro, f3o = i(u(«)2)ro, 70 = |7i - ^(u(°)^)ro, (5o = |(u(")2)ro, «o - §(1 - %ii)(u(o)2)ro 
and g = V(u(")^)/(u(°)^). fl* is assumed to be parahel to the z-axis, ft* = ri*e^ with Q* > 0. 

— . . . . — P — T 

We spht B into its poloidal and toroidal part, B and B , and represent them in the form 

B^ = V X (Ae^) , B^ = Be^ , (129) 

where A and B are two functions of r, z and t, and e^p is the unit vector in (/3-direction. We assume that the differential 

— T — p 

rotation is so strong that concerning the generation of B from B all contributions to S except the /3o-term can be 
neglected. Considering rj again as constant and starting from the induction equation (^ for B we arrive at 

rj^A'A - aon*g,B + ^g ■ V{rA) - {60 - Ko)n*^ " ^ = « 
r oz at 

r^^A'B ~ { — — ^ ——){rA) - — = 0, 130 
^ or oz oz Or' ot 



where 77^ = 77 + /3o and 



A7 = ^(i^(./)) + |^ 

or^rOr ' az^ 



Let us restrict ourselves to a local analysis of these equations, that is, to an investigation in some finite region of 
the rz-plane only. For this purpose we assume that there the coefficients ap, 70, 5o and Kq as well as g are constant. 
Further we assume that there w depends linearily on r and is independent of z, that is, the quantity G = rduj/dr is 
also constant, and duj/dz = 0. We use the ansatz 

{A,B) = Re{(Ao,So)Ji(fc^r)exp(ifc^z + At)} , (131) 

where Aq and Bq are complex constants, Ji is the first-order Besse l function of first kind, kr and kz are real constants 
and A is a complex constant, li gr = this ansatz reduces ( |130| ) to a system of two linear homogeneous algebraic 
equations for Aq and Bq. We are interested in non-trivial solutions only and have therefore to require that the 
determinant of this system vanishes. This leads to 



A = -r^^^k'+i^^ ± \ ~Gn*{iaog.kz + (5o - ^Jo)fc?) - {^^Y (132) 



whith k"^ ~ k1 + k1. If g,. ^ this reduction works only in the limit kr 0, and then (132) applies with fc^ = k^ 



Solutions of the equations (130) with a A possessing a non-negative real part correspond to non-decaying mean 
magnetic fields. We first consider the case gr 7^ 0, in which we have an a-effect. As can be easily seen from ( |132| ) values 
of A with non-negative real part and non- vanishing imaginary part occur if k is sufficiently small. They correspond to 
undamped dynamo waves travelling parallel to the z-axis. Let us further proceed to the case gr = 0, in which there is 
no a-effect. If then —GVI*{5q — kq) > 0, the quantity A is real, and it takes non-negative values for sufficiently small 
k. That is, even in the absence of an a-effect dynamo action proves be possible due to combination of differential 
rotation with the 5 or K-effect. This includes the possibility of such dynamos with J-effect mentioned above. We 
see now that 5 and K-effect are in competition and that they compensate each other if ^0 = ^0- We note that under 
our assumptions (5o — kq = ^(q — l)(u'^"^^). Since g > 1 we may conclude that the dynamo requires G < 0, that is 
dw/dr < 0. 

Of course, the results of the local analysis of the dynamo equations should be confirmed by solving them in all 
conducting space using proper boundary conditions. This has been done so far in the investigations referred to above 
for dynamo models involving differential rotation and ^-effect. 
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9. CONCLUDING REMARKS 



In this paper we have shown a procedure to calculate the mean electromotive force S for a magnetohydrodynamic 
turbulence. The bounds of its applicability result mainly from the use of a closure assumption for the deviation of 
the turbulence from that for zero magnetic field and zero rotation. As explained above it can only be justified for 
sufficiently small mean magnetic fields and small rotation rates of the fluid. For simplicity we have restricted ourselves 
to the case in which there is, apart from the rotation, no mean motion. Specific results have been derived for the 
limit in which the mean electromotive force is linear in the mean magnetic field and a Kolmogorov-type turbulence. 

Let us compare results obtained in the kinematic approach on the basis of the second-order correlation approxima- 
tion, or first-order smoothing, for the high-conductivity limit with results of our procedure. Take as a simple example 
the case of homogeneous and isotropic turbulence. Then we have £ — — aB — /3V x B for sufficiently weak variations 
of B in space and in time. In the kinematic approach under the conditions mentioned we have 

1 ("^ 

a = - (u(x,i) • (V X u(x,t-T)))dT , (133) 
«J Jo 

P ^ II (u(x,i).u(x,t-r)))dT. (134) 

This is often expressed in the form 

a=i(u.(Vxu))ri"),, P = ^{u')r(^l. , (135) 

with properly defined correlation times riorr and riorr- The validity of these results can only be readily justified under 
the condition ut/1 <C 1, where u, I, and r are typical values of the velocity and of the length and time scales of the 
velocity field. This condition, however, is problematic in view of applications to realistic situations, and the validity of 
these results beyond this condition is questionable. Basically it is possible to improve the approximation by including 
higher-order terms in u but this is very tedious. 



The results ( |135| ) have the structure of our result given by ( |116| ), specified by b^"^ = 0, and our result (117). The 
validity of our results, however, is not restricted by a condition like ut/1 1, but only by the applicability of the 
T-approximation (|37|). That is, wc have in any case a much wider range of validity. 

We point out that despite the formal similarity of the mentioned results gained in the kinematic approach and those 



derived here, there is a basic difference between them. In the first case we have origi nally, that is in (|133|) and (134) 



correlations between values of u taken at different times, which are often as in ( 135 ) expressed by u at a given time 
and a correlation time, but in the second case we consider from the very beginning only correlations between values 
of u at the same time. In that sense there is no simple connection between the two kind of results. 

We may use the framework explained above also beyond the limit of very small mean magnetic fields and study, for 
example, a or ^-quenching at least for not too strong fields. Then, of course, Vij and can no longer be replaced 

by v^^^ and mf^'^ . Instead we have to derive equations for Vij and m.^ corresponding to equation (|3C| ) for Xij i insert 
their solutions depending on B in /y in ( ^5|) and follow the above pattern of the determination of Uij and bijk, or a, 
f3, . . .. By the way, then even in the case b^°) = there are contributions of to these coefficients which, of course, 
vanish like ml^^ with B. 



In this context it suggests itself to study in addition to the mean electromotive force £ also the mean pondcromotive 
force !F. It has a part independent of B, which can be calculated with the help of u^-^^ and TOj-^"* only. Its general 
form can be derived on the basis of solutions of the equations for vij and niij mentioned. 
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APPENDIX A: RELATIONS WITH £ijk 



For the derivation of ( p^ it is useful to know the relation 



£ijk''ik + [£iklKj — SjklKi) ^2 ~ ^ijkKk ' K^'') 
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which apphes to arbitrary vectors k and $7. 

In view of the derivation of ( |67| ) we recaU the identity 

— SinSjmSkl — SilSjnSkm — SimSjlSkn ■ (A2) 



APPENDIX B: DERIVATION OF EQUATION (|3o|) 



In the calculations of dxij / dt on the basis of the equations ( |15[) and (|l^) contributions to this quantity occur which 
have, e.g., the form of 



Xij(k,R) = J (5i(u,B;k + K/2)uj(-k + K/2))exp(iK-R) (fK 
i I {kk + Kk/2){ui{k + K/2~Q)uj{-k + K/2))'^k{Q)exp{iK-R)d^K(PQ . 



(Bl) 

In the last expression we may change the sequence of integration so that the inner integral is over K and the outer 
integral is over Q. In the inner integral we may further change the integration variable K into K — Q, denoted by K' 

in the following. In this way, and using QkBk = 0, we obtain 

Xy(k,R) = ^ J {kk + KU2){u,{k- Q/2 + K' /2)uj{-k+ Q/2 

+K'/2))'Bk{Q) exp[i(K' ■ R + Q • R)] d^K' d^Q . 
Remembering the definition of Vij (k, R) we can rewrite this into 

1 fdv,j{k-Q/2,-R)' 



(B2) 



X,,(k,R) 



ifcfcWy(k-Q/2,R) 



dRk 



Bfc(Q)exp (iQ-R) d^Q 



(B3) 



The fact that B varies only on large scales, that is, B is only non-zero for certain small |Q|, suggests to use the Taylor 
expansion 



Wy(k-Q/2,R)-t;y(k,R)- 



1 f dv,,{k, R) 

2 ^ dkk 



This yields 



Xy(k,R) ~ [i{k-B) + ^(B • V)]i;„ (k,R) - fcfei-„,(k,R)Sfc,, , 



(B4) 



(B5) 



where Viji = ■^dvij /dki. According to our assumption a term of the second-order in V was neglected. The contributions 
to the Taylor expansion indicated by O(Q^) only lead to terms of higher order in V and need not to be considered. 
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